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ABSTRACT 

t> : 

We study the influence of general lens and source velocities on the gravita- 
j> ■ tional deflection of light by single and two-point-mass microlenses with general 

axis orientation. We demonstrate that in all cases the lens equation preserves 
its form exactly. However, its parameters - the Einstein radius and the binary- 
lens separation - are influenced by the lens velocity. In Galactic microlensing 
settings the velocity mainly affects the inferred separation for wide binary-star 
or star+planet microlenses oriented close to the line of sight. We briefly discuss 
the case of lenses moving with highly relativistic velocities. 

o : 

£3 . Subject headings: gravitational lensing — relativity 

c§ : 
> : 

"X ■ 1. INTRODUCTION 

The time-dependence of Galactic gravitational microlensing events (Paczyhski 1996) is 
given by the relative motion of the source of light, the lens, and the observer near per- 
fect alignment. These events have been generally successfully analyzed using a quasi-static 
approach, assuming that light is deflected at each instant in a static source-lens-observer 
configuration. This approach is well justified - the lens and source velocities relative to the 
observer are of the order of a couple hundred kms -1 , i.e., ~ 10~ 3 of the speed of light c. 

Nevertheless, with the introduction of image subtraction techniques (Alard & Lupton 
1998) the accuracy of measured microlensing light curves has increased in ideal cases to 
sub-percent levels. In addition there are prospects of high-precision observations of the 
astrometric microlensing effect (see Boden et al. 1998; Han 2001 for the single- and binary- 
lens cases, respectively) by space-based interferometers such as the Space Interferometry 
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Mission (SIM) 1 . Given these developments, one can expect that corrections to the light 
curve or the angular image geometry down to the order ~ 10~ 3 might be detectable. It is 
therefore interesting to avoid the quasi-static approach and investigate the effect of general 
lens and source velocities on light deflection. 

Previous theoretical research has mostly concentrated on the effect of a relative lens to 
observer velocity on the deflection angle by a single lens. Pyne & Birkinshaw (1993) derived 
the first-order effect for a velocity of arbitrary direction, showing that the deflection angle 
increased for lenses moving away from the observer and decreased for those moving towards 
the observer. The same result was later confirmed by Frittelli (2003). Recently Wucknitz & 
Sperhake (2004) presented results for an arbitrarily large but purely radial lens velocity. The 
most general theoretical result can be found in the detailed treatise by Kopeikin & Schafer 
(1999), which includes a derivation of the deflection angle for an ensemble of lenses of general 
velocities. 

In this paper we study light deflection for general lens and source velocities (arbitrarily 
oriented and arbitrarily large) specifically for single and two-point-mass lens microlensing 
events. In §2 we derive the general lens equation and demonstrate the velocity dependence 
of its parameters. In §3 we explore the magnitude of the velocity effects for low velocities 
and illustrate the highly relativistic limit. We conclude in §4 by discussing observational 
aspects and the appropriateness of some of the assumptions. 

2. DERIVATION OF LENS EQUATION 

Our approach is similar to the one used by Klioner (2003) for a single lens in motion. 
We utilize the knowledge of light-deflection formulae in the rest-frame of the lens and using 
Lorentz transformations connect the solution to the rest-frame of the observer. We limit our 
accuracy to the usual first order in deflection angle. For example light rays passing close to 
the components of a binary pulsar are thus beyond the scope of this paper. 

In the case of two-point-mass lenses (hereafter "binary lenses" for brevity - includes 
binary-star and star + planet lens systems) we concentrate on the effect of their center-of- 
mass velocity and neglect the effect of orbital velocity. This is justifiable for sufficiently wide 
binaries (with semi-major axis > 1AU). We return to the case of closer binaries in §4. 

For the purposes of the following calculation we set up the observer rest-frame coordi- 
nates with the origin at the center of mass of the lens at time t' = 0, the z' axis pointing 
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towards the observer, the x' axis in the plane of the sky along the projected binary-lens axis, 
and the y' axis perpendicular to it in the plane of the sky. We denote the distance between 
the observer and the (center of mass of the) lens Dl, the distance between the observer and 
the source plane D s , and the distance between the lens and source planes D LS = D s — D L . 

For simplicity we scale all velocities to the speed of light. We denote the lens and source 
velocities measured in the rest-frame of the observer V and W, respectively. At time t' in 
the described observer rest-frame the observer, the center of mass of the lens, and the source 
are located at 

r' Q = (OA D L ) 

r' L {t') = Vt (1) 
r' s (t') = r' S0 + W(t' + D LS ), 

respectively. In the last expression r' so = (PiD s , #2 As, — D L s) is the source position at 
t' = —D LS . The two-dimensional angle (3 is the angular position of the source in the plane 
of the sky. Note that an undeflected photon arriving at the observer at t' — Dl passed the 
lens at t' — and the source at t' = —Dls- The lens and source velocities and distances are 
thus measured at these retarded times. 

We denote the total mass of the lens M; in the binary-lens case the two lensing bodies 
have masses Ma = H-aM and Mb = /AbM, respectively. As hinted earlier, we place no 
restrictions on the mass ratio of the two lenses - our results thus hold for single- and binary- 
star lenses, as well as for star+planet lenses. We define the coordinates in the rest-frame 
lens of the lens with the origin at the center of mass of the lens, the x axis along the 
physical binary-lens axis, the y axis in the plane of the sky parallel to y', and the z axis 
perpendicular to both. In these coordinates the lenses are located at r A = (— hbR, 0, 0) and 
r B = (haR, 0, 0), where R is the distance between the two lenses (the intrinsic separation). 
We denote the angle between the binary-lens axis and the plane of the sky (, oriented so 
that a small positive value brings lens A closer to the observer. The directions of the x and 
z axis thus coincide with those of the observer's x' and z' only if the binary-lens axis lies in 
the plane of the sky, i.e., ( = 0. 

Denoting the future asymptotic trajectory of a photon (at the observer) in lens-frame 
coordinates 

r po (t) =l + n t, (2) 
where n is a unit vector, the past asymptotic trajectory of the photon (at the source) is 

r ps (t) = r po (t) + Mi ^—^——^ [l-n- n -(l - n) n ] , (3) 

as demonstrated for example in Will (1981) or Brumberg (1991). To obtain the values of the 
constant vectors I and n , we transform the future asymptotic trajectory to observer-frame 
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where the rotation matrix 



coordinates, in which the light ray arriving at the observer is 

r; o (t') = r' o + n' (t'-D L ). (4) 

The unit vector n' — (1 + 6 2 Y l l 2 (— 9± : — 9 2 , 1) describes the direction of the light ray at 
arrival. The two-dimensional angle 6 denotes the angular position of the image in the plane 
of the sky. 

The conversion between the two coordinate systems is given by 

i) ) = MV) ( R(C)U ) ■ (5) 

(cos( sin( \ 
10 (6) 
— sin £ cos £ / 

corrects for the orientation of the binary- lens axis, and the Lorentz boost matrix (e.g., Misner 
et al. 1973) is 

A°'o = 7, A ', = A*'o = 7Vi, A* - = ^ + ^7 - 1) , (7) 

with 7 = (1 — V 2 ) -1 / 2 . From equation (5) we express r po (t) and by comparison with equa- 
tion (2) we get the two vectors I and n describing the future asymptotic light ray in the 
lens frame. 

In a similar way we can take the expression for the past asymptotic trajectory of the 
photon from equation (3) and extend it to the source as follows: 

4))^ a(v) (r(c)U))- (8) 

Note that r' s (t' e ) depends on the source velocity W - see equation (1). To get the photon 
position r ps (t e ) we substitute the previously obtained light-ray vectors I and n into equa- 
tion (3). We first use the time component of equation (8) to convert between the emission 
times t e and t' e in the two coordinate systems, then we use the z component to eliminate the 
emission time altogether. 

From the remaining two equations we can finally express the light-deflection angle 
a = (6 — (3)D S /D LS measured by the observer. In the single lens case we get to first order 
in deflection angle 

o(9) = ^£&z^|. (9) 

c 2 D L V 1 - V 2 v 2 
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This expression is in agreement with the previous results of Kopeikin & Schafer (1999), and 
extends the results of Pyne & Birkinshaw (1993), Frittelli (2003), and Wucknitz & Sperhake 
(2004). In the binary-lens case we get to first order in deflection angle 



a(0) = 



AGM{1 - V z 



c 2 D L Vl - V 2 



Ha- 



6~6j 

e-o A 



0-0b 
0-0 b \ 



(10) 



where 6 a = —^b6ab and 6b = ^a6 ab are the apparent angular positions of the two lenses, 
and the apparent angular separation vector (from lens A to B) in the sky is 



6ab(V) 



1 

~D~l 



l-yi-v 2 

V 2 



VR 



V, 



v~ 



(11) 



Here the tangential lens velocity V± = (V x , V y ) and the vector R = (R±, R z ) = R(()(r B —r A ) 
is the separation vector of the two lenses in their center-of-mass rest-frame at time t' = 0. 
Its component R± lies in the plane of the sky and R z is oriented along the line of sight to 
the observer (for a positive R z lens B lies closer to the observer). 

The first obvious result is the independence of the deflection angle on the source velocity 
W. We note that we have demonstrated this only to the first order in deflection angle, the 
source velocity might have an influence at a higher accuracy. 

The structure of equations (9) and (10) allows us to readily define the velocity- dependent 
angular Einstein radius 



0e(V) 



UGMD ls l-V z 



c 2 D L D s ^/T^V 2 " 



(12) 



We can re-scale all angular quantities by 9e(V) and thus convert from ((3 : 6,6a,6b,6ab) 
to (y, x, x A ,x B , d). The obtained lens equation for a single lens is 



y 



The lens equation for the binary lens is 



x 

X = 2 



y - x = -fi A - 



X 



\X — Xa 



Xa X — Xb 
2 - A*B, 



\x 



XB 



(13) 



(14) 



This interesting result shows that in either case the lens equations have exactly the 
same form as in the static case. When analyzing microlensing observations one can thus use 
exactly the same formulae as in the usual quasi-static approach. However, when interpreting 
the fitted parameters one has to realize that the Einstein radius, and in the binary-lens case 
also the lens positions and their separation, depend on the lens velocity. Instead of #e(0) 
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one obtains 9e(V) and instead of the angular lens separation vector d(0) = R±/[D L 9 E (0)] 
one obtains 

= ^ |d(0) + ( 4 - ^5 4 + y r d(0) ] j A J , (15) 

where d z = R z /[Dl9e{0)} is the ^-component of the binary-lens separation vector (the 
"depth" of the binary lens along the line of sight) re-scaled by the Einstein radius of the 
lens. In the following section we study the velocity effect on these parameters, in order to 
assess its influence on the inferred physical parameters of the microlenses. 



3. VELOCITY EFFECTS ON INFERRED LENSING PARAMETERS 

Light-curve analysis of simple microlensing events does not directly yield the angular 
scale of the event geometry. Event parameters such as the Einstein radius crossing-time and 
the impact parameter are obtained scaled to 9e, and thus will be affected in the same way 
as the Einstein radius. In addition, the scaled lens separation in the binary-lens case will be 
affected as shown by equation (15). 

In events beyond the simple model, such as caustic-crossing (source-transit) events or 
parallax events as well as in events observed astrometrically, it is also possible to measure 
the angular scale. Therefore, in such events the effect on the angular Einstein radius and 
the angular lens separation given by expression (11) is potentially of interest. 

In the following subsections we study the effects of the lens velocity on the angular 
Einstein radius 9e(V), the angular lens separation 0ab{V), and the Einstein-radius scaled 
lens separation d(V). In §3.1 we explore the low-velocity case, which is astrophysically 
relevant for typical Galactic microlensing settings. In §3.2 we demonstrate the results for 
the high-velocity regime. 



3.1. Low Velocities 



The second-order V expansions of expressions (12), (11), and (15) are 



9 E (V) ~ 9 E (0) 



V z 2V? + K 2 



1 -y + 



o AB (v) ~ e AB (o) + i- {r z v ± + 1 -[v z r z - V ± -R ± ] v ± 



d(V) ~ d(0) 



V z V z 2 - 2V? 

1 + — + — - 

2 8 



+ VI 



[1 + V z )d z --V ± .d(0) 



(16) 
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where #ab(0) = R±/D L . The angular Einstein radius has a linear order effect only if the 
lens has a non-zero radial velocity, while any purely tangential velocity produces a non-zero 
second order effect. As expected from previous results, the radius increases for lenses moving 
away from the observer and decreases for those moving towards the observer. However, to 
get a 10~ 3 effect on Oe would require an unlikely radial velocity cV z = 600 kms^ 1 . We 
conclude that the effect of lens velocity on the angular Einstein radius is not observationally 
significant. 

From equation (11) we can see that any purely radial velocity of the lens has no effect 
on the angular separation of the lenses 9ab- The angular separation vector has a linear 
order effect only if the binary lens has a non-zero tangential velocity and at the same time 
its axis is tilted from the plane of the sky. The effect is oriented along the vector of the 
tangential velocity, its direction depends on the sign of R z . The magnitude of the effect 
[0ab(V)/9ab(0) — 1] is R Z (V±-R±)/R 2 ± with a maximum value of V±ta,n( for a purely 
tangential velocity parallel to the projected binary-lens axis. A low value of V± can thus 
be offset by a high angle of inclination of the binary-lens axis. For a tangential velocity 
cV± = 200 kms^ 1 we get a 1% effect for ( = 86°, which corresponds for example to a 
projected separation « 1AU for a binary lens with a physical separation R « 15AU. 
The geometric conditions required for even higher effects are no less realistic. We see that 
the effect is observationally significant for binary-star or star+planet lenses aligned nearly 
along the line of sight. 

To investigate the effect of velocity on the size of the Einstein-radius scaled lens sepa- 
ration vector d(V) we have to treat separately binary lenses oriented along the line of sight, 
for which d(0) = 0. In this case we get for the first-order absolute effect 



For a given total velocity V the maximum effect V\d z \ occurs for a purely tangential velocity. 
The zero minimum effect occurs for a purely radial lens velocity. 

For a general binary lens not oriented along the line of sight d(0) > and the first-order 
absolute effect is 



d(V) - d(0) ~ | d z V ± 



(17) 



d(y)-d(o)^<i(o) + ^-ii4. (is) 




A straightforward computation shows that the maximum value of V \fd 2 {fd) + 4df/2 is 
achieved for a velocity orientation 




(19) 
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The minimum value has the same effect with a negative sign and it occurs for an opposite 
velocity, while a zero effect occurs for velocities perpendicular to velocity (19). It is interesting 
to note that the maximum effect on the scaled lens separation occurs generally for a different 
velocity orientation than the maximum effect on the angular lens separation. Figure 1 
illustrates the dependence of the absolute first-order effect on the direction of the lens velocity 
for different orientations of the binary-lens axis. In terms of lensing parameters the maximum 
first-order absolute effect for a binary lens with its axis tilted by ( from the plane of the sky 
is 

i(V) _ , (0) * o. 01 ( cV \ /Wt 1 (R\ (D^\ 1 

v ; v; 2 \200kms- 1 ) \lmas) \60AUJ\4kpcJ v ; 

for a velocity direction given by equation (19). The effect for higher lens velocities, wider 
intrinsic lens separations R, and lower Einstein radii can thus be a significant fraction of the 
Einstein radius. The maximum relative effect [d(V)/d(0) - 1 ] is V ^cos- 2 C - 0.75, which 
is largest for ( — > 90°. In this regime it coincides with the maximum relative effect for the 
angular separation Q ab derived above. 

To summarize the low- velocity results, first-order effects are caused by the radial velocity 
of the lens and/or the "depth" of the binary lens along the line of sight. In the case of a 
single lens the effect of lens velocity on the inferred lensing parameters is not significant. In 
the case of binary-star or star+planet lenses the inferred values of the angular and Einstein- 
radius-scaled lens separations (as well as the linearly related lens positions) can differ from 
the values at zero velocity by more than 1% if the axis of the lens system is oriented close to 
the line of sight. In particular, the effect on the scaled lens separation can be a significant 
fraction of the Einstein radius mainly for lenses with wide intrinsic separations. 

3.2. High Velocities 

The general dependence of the angular Einstein radius on the lens velocity as given by 
equation (12) is illustrated in Figure 2. We can see that the high- velocity regime depends 
on the orientation - the case when the lens moves directly towards the observer has to be 
treated separately. 

If we increase the velocity V — > 1 in any other direction (including directly away from 
the observer), the angular Einstein radius eventually diverges 0e(V) oo. For the angular 
lens separation equation (11) gives us a finite result, and the scaled lens separation thus 
vanishes d(V) — > 0. Hence, such a "binary lens" would in effect behave like a single lens. 

We note that in this regime the approximation of first order in deflection used in this 
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paper eventually breaks down. However, the divergence of the Einstein radius is very slow - 
on the order of ~ (1 — V)^ 1 ^. From Figure 2 we can see that even for V = 0.9 the Einstein 
radius is at most greater by a factor of two (if the lens moves directly away). Even with a 
factor of 100 or more the approximation remains valid. It is thus sufficient to formally ignore 
the results for V = 1. 

A lens moving with a velocity V — > 1 directly towards the observer will have a vanishing 
angular Einstein radius 9e{V) — > 0. The angular lens separation 0ab — * R±l Dl - the lenses 
appear to be at their "true" positions. The scaled lens separation thus diverges d(V) — > oo. 
The decreasing Einstein radius means that such a "binary lens" would behave like two single 
lenses with decreasing strength, as the lens moves together with the arriving photons. 

4. DISCUSSION 

In §3.1 we have demonstrated the low- velocity effects for different binary-lens orienta- 
tions and intrinsic separations. Nevertheless, from the observational perspective we must 
take into account the fact that not all binary-lens configurations can lead to microlensing 
events detectable as two-point-mass lens events. While this depends on the exact source 
trajectory, from a statistical point of view the Einstein-radius scaled lens separation d plays 
the main role. If the separation is too small or too large, most events would appear as 
single-lens events. 

For the event to be detectable as a two-point-mass lens event the lens separation has to 
fulfil d(0) G (d mini d max ) . This leads to limits on the binary-lens axis orientation 

cosCe(rfm „,,^. )x (^)- 1 (_^)(M21). (21) 

The values of the binary separations from the 21 binary events published by the MACHO 
team (Alcock et al. 2000) range from 0.421 to 2.077 with an outlier value of 7.454 (a possible 
binary-source event). The 18 events detected by OGLE-II (Jaroszyhski 2002) have values 
from 0.355 to 2.917, and the values for the 15 events detected by OGLE-III (Jaroszyhski et 
al. 2004) range from 0.352 to 3.457. 

To obtain a rough estimate we set d m i n = 0.3 and d max = 4. Keeping Di and &e fixed at 
the values used in equation (21), for an intrinsic binary separation R = 15AU we get limits 
on the axis angle C G (0°, 85°), for R = 60 AU we get C e (75°, 89°), and for R = 240 AU we 
get C G (86°, 89.7°). In the first case we have a wide range of possible orientations - however, 
if we approach 90° where the velocity effects are strongest, it would be difficult to detect such 
an event as a binary- lens event. The two cases with higher intrinsic separations demonstrate 
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that such lenses can have stronger velocity effects while being detectable as binary events, 
albeit for a narrower range of axis orientations. 

In this work we concentrated on the center-of-mass velocity of the binary lens and 
neglected its orbital velocity. However, for binaries or star+planet systems with semi-major 
axis <C 1AU the assumption of small orbital velocity breaks down. The results of this paper 
indicate that the velocity effects are proportional to the intrinsic lens separation, being 
caused by the lens motion during the passage of the light ray in its vicinity. While the 
orbital velocity grows as R~ x l 2 with decreasing R, its product with the separation decreases 
as R 1 / 2 . We conclude that the effects of orbital lens velocity on light deflection are not 
significant even in binary-lens systems closer than 1AU. 

This work was supported by grant GACR 205/04/P256 from the Czech Science Foun- 
dation. 
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z $ = 90° 



Fig. 1. — Radial plots of velocity effect on binary-lens separation d(V) —d(0) given by equa- 
tions (17) - (18), depending on velocity direction in the observer+binary plane (z axis points 
to observer; x axis in the plane of the sky). Individual plots are for different binary-axis ori- 
entations C (dot-dashed lines: binary axis). Effect radially scaled in units of V \/d 2 (0) + d 2 z ; 
arrows mark direction of maximum positive effect given by equation (19). Three-dimensional 
extensions of the plots are rotationally symmetric: for ( ^ 90° around the arrow, for ( = 90° 
around the z axis. 
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Fig. 2. — Contour plot of angular Einstein radius ratio 9e(V)/9e(0) as a function of radial 
and tangential lens velocities V z and V± (V z positive for motion towards observer). Plotted 
contour values from 0.5 (lower right) to 2 spaced by 0.1: dashed for values < 1; solid bold 
for 1; solid for values > 1. Dot-dashed curve: V = 1 (speed of light). 



